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We propose an analytical model devoted to explain the anisotropy of the electrical resistance ob- 
served below the critical temperature in granular high-T c superconductors submitted to a magnetic 
field H. Reported experimental results obtained on a YBCO sample show that the superconducting 
transition occurs in two stages, with a steep drop of the resistance at T c and a subsequent, smoother 
decrease. In this second stage, the resistance vs. temperature curve is strongly dependent not only 
on the field intensity, but also on the angle between H and the macroscopic current density j. 

We start from the assumption that the resistance below T c is mainly due to the weak links between 
grains. In the model, weak links are thought of as flat surface elements separating adjacent grains. 
We calculate the probability for a weak link to undergo the transition to the resistive state, as 
a function of the angles it makes with the external magnetic field H and the macroscopic current 
density j. In doing this, an important role is given to the strong non-uniformity of the local magnetic 
field within the specimen, due to the effect of the screening supercurrents flowing on the surface of 
the grains. Finally, we calculate the electrical resistance of the sample in the two cases H _L j and 
H || j. The predictions of this simple model turn out to be in reasonable agreement with reported 
experimental results obtained on a YBCO granular specimen. 

PACS numbers: 74.25.Fy; 74.50.+r 



I. INTRODUCTION 



If a granular sample of a high-T c superconductor 
(HTSC) is cooled down to its critical temperature (let 
us call it T c o) in the presence of a magnetic field, its elec- 
trical resistance suddenly falls to a value that can be as 
low as 30% of the normal-state resistance just above T c o 
(see Fig. 5 in Section IV). This effect, common to many 



cuprates, is due to the superconducting transition of the 
grains while the iiatergrain regions (weak links) remain in 
the resistive statcfil. On further cooling, the sample resis- 
tance gradually decreases, eventually becoming zero at a 
temperature T c whose value depends on the applied mag- 
netic field. In this second stage of the superconducting 
transition, the transport properties of the specimen are 
entirely controlled by the weak links between grains, that 
can be thought of as S-N-S Josephson junctionscl with 
randomly-distributed critical energies. According to this 
approach, when the temperature is lowered the super- 
conducting wavefunctions of the grains gradually lock in 
phase. This gives rise to long-range coherence and finally 
to the bulk superconductivity. Of course, the transition 
of each junction from the superconducting to the resis- 
tive state (or vice-versa) is controlled by the tempera- 
ture, the current density crossing the junction and the 
local magnetic field. These three parameters, together 
with the distribution of the Josephson critical energies, 
are thus expected to play a major role in determining 
the electrical behaviour of the material. The transition 
also depends on wether the magnetic field is applied dur- 
ing the cooling of the sample (field cooling, FC) or after 



cooling (zero field cooling, ZFC) . Actually, in the FC case 
field penetration and trapping within the grains reduces 
the transition probability of the weak links, and conse- 
quently the value of the electrical resistance, by reducing 
the flux compression factor. This last quantity will be 
introduced and discussed in Section III of this paper. 

Finally, in many cases, the electrical properties of gran- 
ular HTSC below T c q are found to depend on the angle 
between the magnetic field H and the macroscopic cur- 
rent density j, even if the grains are randomly oriented. 
This indicates the existence of an electrical anisotropy of 
the material, induced by the application of an external 
magnetic field. That also this anisotropy is related to 
the complex transition dynamics of the weak link net- 
work is demonstrated by its almost complete absence in 
high-density polycrystalline MgB2 samples, where grains 
are connected through metallic contacts and there are no 
weak linksa. 

In the following, we will propose an explanation for the 
observed electrical anisotropy in a zero-field-cooled gran- 
ular HTSC. The leading idea of the model is that the 
screening supercurrents flowing on the surface of the su- 
perconducting grains create a local "demagnetising field" 
that adds to the external one creating a strongly non- 
uniform field distribution in the integrain regions. We 
will show that this field distribution gives rise to a struc- 
tural anisotropy in the network of resistive weak links, 
so that the material behaves as a uniaxially anisotropic 
medium for the current transport. By starting from sim- 
ple hypotheses, which will be discussed and supported by 
experimental evidences, we will develop a simple network 
model that allows calculating the electrical resistance of 
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the material as a function of the magnetic field and of the 
current density in the two cases where H _L j and H || j. 
Finally, we will compare the results of our calculations 
with an extensive set of experimental results obtained on 
a YBCO granular specimen. 



II. THE APPROACH TO THE PROBLEM 

Let us focus on a zero-field-cooled (ZFC) granular su- 
perconductor. Suppose to feed it with a current of den- 
sity j and to apply a magnetic field H such that the bulk 
superconductivity is disrupted, but the grains remain in 
the Meissner state. As long as the grains exclude the 
magnetic field, experiments show that the resistance vs. 
magnetic field curves are very nearly reversible. This re- 
versibility shows that it is the flux pinning within the 
grains that originates hysteresis, while flux pinning in 
the intergrain regions is negligible. In these conditions, 
the electrical resistance is found to depend on whether 
the magnetic field H and the macroscopic current den- 
sity j are parallel or perpendicular to each other. This 
anisotropy was already evidenced by some measurements 
of critical GurrentoQ, resistivityffi], magnetizationEl, power 
dissipationtl and I-V characteristics!^ in different HTSCs. 
Most of the relevant papers propose an explanation for 
the anisotropy based on the conventional theories of the 
current-driven vortex motion in the mixed state. Within 
this picture, the anisotropy arises from the fact that the 
Lorentz force between transport current and vortices de- 
pends on the angle between H and j. Actually, this re- 
quires that the material behaves as a nearly homogeneous 
medium with an effective penetration depth A e ff, where 
the magnetic flux penetrates in the form of vortices as 
in conventional type-II superconductors. As pointed out 
by Ginzburg et alt3 this approach is reasonable as long 
as A e ff is much greater than the average grain size, and 
fluxons can be taken parallel to the external field. Even 
in this case, however, the conventional theories must be 
improved - for example to explain why, in the "H || j" 
case, the voltage drop across the specimen is far from 
being nearly zero. Finally, experimental results that will 
be presented in the following section clearly indicate that 
the anisotropy gradually vanishes when the temperature 
approaches T c o , and this behavior is not easily explicable 
within the picture described so far. 

In the following we will propose a completely differ- 
ent explanation for the electrical anisotropy induced in 
granular HTSCs by the magnetic field. We will start by 
representing the material as a set of irregularly shaped 
grains connected through thin intergrain regions that be- 
have as weak links. We will think of these weak links 
as resistively shunted Josephson junctions with a per- 
fect ohmic behavior in the normal state - that is, with 
a normal-state resistance independent of both magnetic 
field and current, as well as of the angle between them. 
This assumption finds again support in the fact that, at 
temperatures close to T c q - when almost all the weak 



links are in the resistive state - the anisotropy disap- 
pears. This suggests that the anisotropy is related to the 
spatial distribution of superconductive and resistive weak 
links within the material, rather than to an intrinsic de- 
pendence of the intergrain resistivity on the orientation 
of H and j . 

The problem now is to understand why the distribu- 
tion of resistive weak links should be anisotropic. We 
argue that this distribution is mainly determined by the 
distribution of the local magnetic field intensity in the 
intergrain region. As a matter of fact, the local magnetic 
field Hi is given by the superposition of the external field 
H and of the magnetic field created by the screening su- 
percurrents that flow on the grains' surfaces. As a result, 
the local field intensity Hp in a given weak link can be 
very different from H - and can vary very drastically 
from a weak link to another. Simple geometrical con- 
siderations lead to the conclusion that this variation is 
related to the spatial orientation of the weak links with 
respect to H. The strong non-uniformity of the local 
field makes the transition probability of the weak links be 
anisotropic. Therefore, the spatial distribution of weak 
links that undergo the transition to the resistive state 
becomes anisotropic as well, and the material behaves as 
a uniaxially anisotropic medium for the current. 

The influence of the screening supercurrents of the 
grains on the intergrain region was already invoked 
to explain other interesting properties of granular su- 
percondoictoES | _such as the hysteresis of the critical 
currentli3'l!j'oE3 and the ac magnetization cucwsB. In 
the original approach by Evetts and Glowackill3 it was 
assumed that a sufficiently weak magnetic field can be 
excluded both by grains and by superconducting "is- 
lands" bounded by closed paths, called "rings", made up 
of grains connected through weak links with relatively 
high critical currents. According to their discussion, the 
screening supercurrents flowing along the boundaries of 
the superconducting regions create "flux compression" 
in the surrounding weak links where the field has pene- 
trated, which thus experience a magnetic field more in- 
tense than the applied one. It is worthwhile to notice 
that, in that paper, the local magnetic field outside the 
superconducting regions was supposed to be everywhere 
greater than the applied field. 

On the contrary, experimental studies of the ac magnet 
tization of granular HTSCs led Chandran and ChaddahEl 
to suggest that the screening supercurrents flowing on 
the grain surface give rise to "flux compression" in the 
weak links laying on planes parallel to the external field 
H, and to an almost complete magnetic shielding of the 
weak links perpendicular to H. In spite of the oversim- 
plification implicit in this model - grains are thought to 
be cubic and arranged in a regular lattice, as in Fig.l - 
the idea it is based on can be safely assumed to explain 
the origin of the magnetic field-induced anisotropy of the 
resistance in granular superconductors. 

In the present paper, we shall neglect the possi- 
ble contribution of superconducting "rings" to the non- 
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uniformity of the local magnetic field in the intergrain re- 
gions. Actually, the experimental results to which we will 
compare the theoretical predictions of our model were ob- 
tained on a YBCO granular specimen with a small crit- 
ical current density (less than 10 5 A/m 2 at T — 27 K). 
A simple calculation shows that the maximum magnetic 
field created by this current density flowing on a circular 
ring made up of grains and superconducting weak links 
is definitely negligible with respect to the magnetic fields 
considered here,-ejien if the ring is very small. Notice 
that also in refsEfrtj experimental results were reported, 
suggesting that flux trapping or exclusion is mostly due 
to the grains, rather than to persistent superconducting 
loops in the weak link network. 

The model we are going to present in the following 
arises from a generalization of the idea by Chandran and 
ChaddahQ to a more realistic case, in which the grains 
have irregular shape and the weak links are randomly 
oriented in space. We will show that this model quan- 
titatively explains the results of a set of resistance mea- 
surements we carried out on a YBCO granular specimen. 



III. THE MODEL 



(a) 



(b) 




external magnetic field 
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magnetic field due to 
the screening supercurrents 



(c) 



weak link 
with 9=0 



H 



weak link 
with 9=7t/2 



A. The simplest case 

To discuss the origin of the field-induced resistance 
anisotropy in granular HTSCs, let us start with the anal- 
ysis of the ideal, simplified case in which the material is 
made up of a set of identical cubic grains arranged in 
a regular lattice, as shown in Fig. [l[ Within this sim- 
ple picture, the weak links are represented by the flat, 
square surfaces separating adjacent grains. Let the mag- 
netic field H be applied parallel to one of the grain edges, 
for example along the vertical direction. As previously 
discussed, H is assumed to be intense enough to destroy 
the magnetic screening of the sample as a whole (due 
to the supercurrent flowing on the sample surface), but 
weak enough not to penetrate into the grains. Let 9 be 
the angle between the field H and the normal n to a 
given weak-link surface. It is clear that, in this simple 
model, only the values 9 = (corresponding to n || H) 
and 9 = tt/2 (corresponding to nl H) are possible. 

The screening supercurrents flowing on the surface of 
each grain, which is supposed to be in the Meissner state, 
create a "demagnetizing field" that cancels out the exter- 
nal magnetic field H inside the grain. In the surrounding- 
weak links, the local magnetic field is given by the su- 
perposition of the demagnetizing fields of adjacent grains 
and of the external field H. As shown in Fig.l, the re- 
sulting local magnetic field is much more intense than the 
external one (i.e., the flux is "compressed") in the weak 
links having 9 = tt/2, while it is zero in the weak links 
having 9 = 0. 

In the hypothesis that the weak links behave as Joseph- 
son junctions, their transition to the resistive state occurs 
when the density of the current crossing the junction is 



FIG. 1: In the simplest model for a poly crystalline HTSC, 
grains are cubic, arranged in a regular lattice and separated 
from one another by intergrain regions that behave as weak 
links, (a) If a weak magnetic field (thick arrows) is applied 
to the zero-field-cooled material, the screening supercurrent 
flowing on the surface of each grain creates a demagnetizing 
field (thin loops with arrows), (b) The effective field in the 
intergrain region results from the superposition of the demag- 
netizing fields of different grains and of the applied magnetic 
field, (c) The resulting field pattern is such that the magnetic 
field is zero in the weak links with 8 = 0, and is greater than 
the applied one in the weak links with 9 = tt/2. This field 
enhancement is usually referred to as "flux compression". 



greater than a critical value j c (T, H) which depends on 
the temperature and on the local magnetic field. 

Therefore, when the current density j is perpendicular 
to the external field H, the condition for the resistive 
transition is easily satisfied in the weak links having 9 = 
7r/2, which are crossed by the current and submitted to 
a strong magnetic field. The resistance of the sample is 
thus different from zero. 

On the contrary, when j is parallel to H the weak links 
having 9 = remain in the superconducting state even 
though they are crossed by the current (provided that the 
current density is not too large) because of the magnetic 
screening of the grains. Therefore, each vertical column 
of interconnected grains behaves as a superconducting 
path for the current, and the specimen resistance drops 
to zero. 

To summarize, the non-uniformity of the local mag- 
netic field makes the spatial distribution of resistive weak 
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links depend on the direction of j with respect to H. The 
macroscopic result is a field-induced anisotropy of the 
transport properties of the material (in particular, of the 
resistivity). 

After this simple explanation of the basic mechanism, 
the generalization of the model to a more realistic situ- 
ation, in which grains have irregular shape and size, has 
now to be considered. We shall adopt a statistical point 
of view and make some simple assumptions, consistent 
with the experimental conditions. 



B. The hypotheses 

Let us identify the weak links with flat elements ap- 
proximating the surface separating adjacent grains, with 
random orientation in space and average area As. Let n 
be the unit-length vector normal to their surface, 9 the 
angle between n and the applied magnetic field H, and f3 
the angle between n and the macroscopic current density 
j . Let Hi and ji be the local magnetic field and cur- 
rent density within a weak link. Notice that, as long as 
the grains are in the Meissner state, must be parallel 
to the weak-link surface - that is, tangent to the grain 
boundary. We will further assume that the weak links 
behave as ideal resistively shunted Josephson junctions, 
with a perfect ohmic behavior above the transition, and 
that they have all the same resistive-state conductance 
per unit surface, g. 

In principle, each weak link undergoes the transition 
from the superconducting to the resistive state when the 
local current density ji is greater than a critical value 
ji^ c (Hi,T). However, complex transient phenomena oc- 
curring at the beginning of the conduction process make 
it very difficult to determine the actual spatial distribu- 
tion of the resistive weak links. These phenomena are due 
to correlation effects acting on the local current distribu- 
tion, which, for instance, prevent the weak links embed- 
ded in a superconducting region to undergo the resistive 
transition. Anyway, when the whole specimen becomes 
macroscopically resistive, a stationary situation similar 
to that sketched in Fig. || must be reached. At the equi- 
librium, the specimen must be thought of as divided in a 
set of equipotential regions separated by resistive layers 
extended throughout the specimen cross-section. 

The equipotential regions consist of several supercon- 
ducting grains and of the weak links between them. Of 
course, all these weak links must be in the superconduct- 
ing state, since any potential drop within these regions is 
forbidden. It is worthwhile to notice that the magnetic 
flux may penetrate within the weak links embedded in an 
equipotential region, but the local current density should 




(a) 



Equipotential 
regions 



(b) 



H 



FIG. 2: Longitudinal cross section of a zero-field-cooled cylin- 
drical poly crystalline HTSC some time after the application of 
a magnetic field H intense enough to make it macroscopically 
resistive. At the equilibrium, the sample consists of equipo- 
tential superconducting regions (gray) , separated by resistive 
layers (irregular solid lines) made up of resistive weak links. 
The effect of the relative orientation of H and j is made clear 
by comparing (a) with (b). It is related to the fact that the 
probability for a weak link to undergo the transition to the 
resistive state is larger when H is parallel to its surface, and 
smaller when H is perpendicular to it. 



remain below its critical value to avoid the resistive tran- 
sition. 

The resistive layers, instead, are made up of resis- 
tive weak links. Since these layers separate two regions 
at different potential, the local current density that 
crosses them must be always perpendicular to their sur- 
face, while the local magnetic field Hf is tangent to it. 
Since the weak links are assumed to behave as ideal 
shunted Josephson junctions with a perfect ohmic be- 
haviors above the transition, the electrical resistance of 
each layer may be considered as independent of both 
and j£. As already pointed out, this last assumption is 
supported by the resistance vs. temperature curves re- 
ported in Fig. 5, measured in a YBCO granular specimen 
described in Section 
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1 Here and in the following we shall use j to indicate the vecto- 
rial average of the current density within the whole specimen. 
Instead, the local current density will be indicated by j^. 



At a given temperature, the transition of a weak link to 
the resistive state is completely determined by Hi and j^. 
As a good approximation, we can say that the transition 
occurs when 



jt = je,o(T) 



Hp 
ttHi 



(1) 
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where je,o(T) is the value of the local critical current 
density in zero magnetic field, and H Q is given by 



H 



W(T)i? g 



Here A(T) is the magnetic penetration depth, R g is the 
mean radius of the grains, and 0o is the flux quantum. 
Equation (|l|) represents the envelope of the Fraunhofer- 
like I vs. H curve of a single Josephson junction in 
the presence of a magnetic field parallel to the junction 
itselftl. Using the envelope instead of the true function is 
usual when a statistical approach is needed - i.e., when 
a large number of junctions enter into the model. Actu- 
ally, as far as the inverse proportionality between je and 
Hi is concerned, the validity of Eq. (|l|) is supported by 
experimental results that will be discussed in section RM. 

For any given value jt of the current density, Eq. (Uf) 
can be interpreted as a condition on the intensity of the 
local magnetic field, He- The critical value of He giving 
rise to the transition will be indicated in the following by 



H c (j e ) = H 



irji 



(2) 



With reference to the equilibrium situation described 
in Fig. 2, the local current density in a resistive weak link 
is given by 



= j • n = jcos/3 



Thus, Eq. (0) becomes: 



H c (je) = H 



irj 



cos (3 



Hc(j) 



cos (3 



(3) 



We will assume that, for a given value of the local cur- 
rent density, the critical fields H c (je) of the weak links fol- 
low a gaussian distribution with mean value < H c (je) > 
and standard deviation a c . In fact, the fluctuation of 
H c {je) around its mean value - within the weak-link en- 
semble characterized by a given value of je - is due to 
several uncorrelated causes: the grain orientation mis- 
match, impurity segregation at the grain boundaries, 
non-stoichiometric local oxygen content, etc. In these 
cases the assumption that the fluctuation is gaussian is 
generally accepted. We further assume that <r c is propor- 
tional to < H c (jg) >. This assumption can be justified by 
observing that when < H c (je) > is, for instance, reduced 
as a consequence of an increment of jj>, also the fluctua- 
tion around its value must change accordingly. Actually, 
it must be noticed that < H c (je) > depends on ji and 
thus on the spherical angles 9 and (p, and that the av- 
erage is intended to be made over the set of weak links 
characterized by given values of these angles. 

Equation (^) indicates that the critical field of each 
weak link depends on the angle (3 between its normal n 
and the current density j. However, the intensity of the 
local magnetic field in the weak link, He, is expected to 
depend in some way on the angle 9 between n and H. 



That this dependence should exist is suggested by the 
simple model of cubic grains, where the value of the lo- 
cal field was easily determined for all the allowed values 
of 9. In the more realistic case of irregular-shaped grains 
we are facing here, flux lines meander through the sam- 
ple, without violating the requirement of continuity. It 
is thus very likely, for instance, that flux lines are forced 
to pass through weak links with a very low value of 9, or 
that weak links with 9 close to ir/2 are almost completely 
screened. Since the value of He in a given weak link char- 
acterized by an angle 9 also depends on the position and 
on the angle distribution of the nearby weak links, it is 
reasonable to assume also in this case a gaussian distri- 
bution of He (6) around its mean value < He(9) > with 
a variance oe proportional, as in the previous case, to 
< He(9) >. Since He should always be tangent to the 
grain surface, we can rather safely assume that the de- 
pendence of < He{9) > on 9 is expressed by the equation: 



< H e (9) >= H'sin9 



(4) 



analogous to the expression of the local magnetic field 
intensity on the surface of a superconducting sphere im- 
mersed in a uniform magnetic fieldEj. In this equation, 
H 1 is a constant magnetic field intensity that is related to 
the external magnetic field through a "flux compression 
factor" k that takes into account the effect of the flux 
exclusion by the grains. Finding an explicit expression 
for k will be the aim of the following section. 



C. Flux compression factor k 

Let's now assume that the magnetic field H is applied 
parallel to the z axis (whose direction is defined by the 
unit-length vector u z ). To determine the value of H', we 
first write down an expression for the mean value of the 
z-component of over the ensemble of weak links with 
the same 9. With reference to Fig. || one finds: 

<Ht <z (6)> = <He(9)sm9> (5) 
= < He{9) > sin9 = H'(sin9) 2 

The mean value of the z component over all the weak 
link ensemble is now given by a simple solid angle aver- 
age: 



< H t ,z > = « H t , z (9) » 
1 

47T 



(6) 



#'(sin( 



2 dn = -H' 



In order to find the compression factor, we express the 
same quantity in a different way. In fact, if no supercon- 
ducting rings survive, the magnetic field is only excluded 
by the grains and its dispersion at the sample edges is ex- 
pected to be negligible. In these conditions, the magnetic 
flux <f> is very nearly constant all along the specimen and 
obviously equal to H-Afj, where H is the intensity of the 
applied magnetic field H (which is taken parallel to the z 
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FIG. 3: Representation of the vectors characterizing the po- 
sition of a single weak link. The figure also shows that the 
local magnetic field is always parallel to the weak-link surface. 
8 is the angle between the unit-length vector normal to the 
surface, n, and the applied magnetic field H. 



the resistive state when the local magnetic field within 
it, Hi(9), becomes equal to a current-dependent critical 
field H c (ji). The transition probability for this weak link 
is thus: 



P tT (9,(3) = f(H t \0)P(Ht>H c (jt))dHt (11) 



axis) and Ah is the specimen cross-section perpendicular 
to it. On the other hand, the explicit calculation of the 
flux gives: 



■dS 



H e , z -dS=< H l>% > A' H (7) 



The integral has been restricted to the portion A' H of 
the specimen cross-section in which the magnetic field has 
penetrated. Provided that there is no flux penetration 
into the grains, A' H can be taken as a constant quantity. 
Finally, < _f/^ z > is obviously the mean value of Hg^ 
over A' H . 

H ■ Ah gives 
(8) 



Comparing Eq.(Q) with the equality $ 

A H 



< H Lz >= H ■ 



A 



H 



and, by comparing this result with Eq. (^) , an expression 
for H' is finally obtained: 



H' 



3 A 



H 



H = kH where k = 



3 A 



H 



2 2 



(9) 



In conclusion, taking into account Eq. ([|), the mean 
value of the local field intensity obtained by averaging 
over all the weak links with the same 9 can be written: 

< H e (9) >= kH sin 9. (10) 



D. Transition probability for the weak links 

As previously pointed out, we suppose that a weak 
link with a given value of 9 undergoes the transition to 



where f(Ht\9) is the distribution of the local field inten- 
sity for a given 9, and P(Hi > H c (ji)) is the probability 
for Hg to be greater than H c {ji). According to previous 
assumptions, f(He\9) can be written: 



f(H e \9) = 



1 



2-KCfl 



exp 



(H e - kH sin ( 



(12) 



where Eq. ( |i0| ) has been used to express the mean value 
of He(9). Similarly, the expression for P(He > H c (je)) 
reads: 



P(Hi > H c (jt)) 



H, 



1 (Hc-<ffc(Ji)>r 

5^ 

e 2<7 c 



2-KCJr 



1 +erf 



dH r 



H f ~ 



<Hcjfy> 

cos f3 



V2a c 



(13) 



where Eq. (|3j) has been used to express < H c (jA > in 
terms of the angle /3. It can be noticed that Eq. ( pT| ) al- 
ready contains all the information about the field-induced 
anisotropy we are dealing with in the present paper, in 
the sense that P tT will have a different expression as the 
current density j is parallel or perpendicular to the ap- 
plied magnetic field H. As a matter of fact, the relation- 
ship between f3 and 9 is different in the two cases. With 
reference to Fig. |], it is clear that: 

if H || j cos/? = cos(9 (14) 

if H _L j cos (3 = sin 9 cos <p , (15) 

where ip is the angle between j and the plane containing 
both H and n. Taking into account Eqs. (14) and Jl5|), 
we shall indicate the transition probability functions for 



7 



H || j and H _L j as Pf (0) and Pf (0,<p), respectively. 
By using Eqs. (|ll]) , (12) and (|l^), and by taking into 
account the expressions for cos [3 given in Eqs. (n4j) and 
(|l5|), one obtains 



Pf(0) = 
and 



27T(Tf 



exp 



{He - kH sin i 



2aj 



1 + erf 



<HcU)> 



V2a c 



dHv 



(16) 



exp 



(i^ - kHsmO) 2 



2o\ 



1 + erf smficosf 



dfl-*. (17) 



E. Calculation of the specimen resistance 

As discussed above, when the specimen resistance is 
different from zero, the specimen itself can be described 
as a series of equipotential stripes, separated by thin lay- 
ers made of resistive weak links. The dependence of the 
transition probability of the weak links on their angular 
position makes these resistive layers look very different in 
the two cases H || j and H _L j, as shown in Fig. |[ Since, 
in agreement with the experimental results reported in 
Fig. g, all the weak links are assumed to have the same 
resistive-state conductance per unit surface g, the con- 
ductance of each layer is simply proportional to its area, 
S. The mean value of the area S can be easily calcu- 
lated by taking into account that each layer is made of 
resistive weak links, and that its projection on a plane 
perpendicular to the current density j must be equal to 
the specimen cross-section area Aj : 



<S>= 



< | COS/3| >re 



(18) 



Here, the subscript "res" means that the average is made 
over the ensemble of the resistive weak links. Thus: 



< | cos/3| > ros = 



j P tr (9,P)\cosf3\dn 



(19) 



where P tr is given by Eq.(|l6|) or Eq.([l7|) according to 
whether H || j or H _L j, and the integrals are extended 
to the whole solid angle Vl — An. The mean value of the 
resistance of each layer is thus 



< S > g 



< | COS/3| > r 



(20) 



Let 71l be the number of resistive layers: the resistance 
of the specimen is thus given by 



R = n L r = n L 



< | cos/3| > r 



(21) 



To evaluate riL, we must first calculate the total num- 
ber of resistive weak links, N, involved in the creation 
of a given pattern of resistive layers. N is proportional 
to the integral of the transition probability of the weak 
links, P tr , over the whole space. More precisely, if Nt is 
the total number of weak links in the whole specimen, N 
is given by: 



N = N T P tr {9,P)dn. 



(22) 



On the other hand, the number of weak links that com- 
pose a single resistive layer, n w , is: 



<S> 



As As- < | cos/3| > rc 



(23) 



where As is the average area of a resistive weak link and 
is of the order of the square of the grain radius R g . Thus 
the number of layers results 



n = — = N — < * C0S/ ^ >rc 

7i W Aj 



and the specimen resistance becomes: 

r As(< | cos 0| > rcs ) 2 



R = N- 



9^ 



(24) 



(25) 



It is worthwhile to notice that, in Eq.(|25|), the speci- 
men resistance turns out to be proportional to the total 
number of resistive weak links. Actually, as is well known 
from the percolation theories applied to the supercon- 
ducting transitiorO, this is not true in the proximity of 
the percolative threshold, since below a minimum num- 
ber of resistive weak links not a single resistive layer is 
generated. This implies that below a given value of H 
the specimen is in the superconducting state and obvi- 
ously Eq. (H|) is not valid. This point will be taken into 
account in comparing the theoretical results with exper- 
imental data. 



F. Anisotropy of the resistivity 

As previously pointed out, all the information about 
the anisotropy of the resistance is already contained in 
the transition probability (see Eq. (jlj) and (p7|)). P tr 
enters directly into the calculation of the specimen resis- 
tance through the mean value of | cos /3| over the ensemble 
of resistive weak links. We shall use R\\ and R± to indi- 
cate the resistance in the case where H || j and H _L j, 
respectively. According to Eqs. (p5|), ( p2] ) and (|l9|), and 
by u sing the expressions for cos/3 reported in Eqs. (|lj) 
and (|15|), one finds: 



As 
1% 



ff (0)dfi ■ 



/ P | | ;r (0)|cos6'|dfi 



(26) 
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and 



R± = 



As 



(27) 



/ F£(0,<p)|sin0cosp|dfi 



fp?(e,tp)dn 



In order to evaluate the anisotropy of the resistance 
(that is, its dependence on the respective orientation of 
H and j) we shall define a parameter r\ such that: 



V 



R± 



£1 
p± 



(28) 



The last equality holds because the current always flows 
in the same direction with respect to the specimen (and 
the direction of the magnetic field is changed instead). 
On account of Eqs.fl26|) and (|27|), 77 can be written as 
follows: 



fp^(e^)dn 

JP^(9)dQ 



JP^(0) I cos 6 1 dn 
J P?(9,tp) I sin 9 cos tp\ dn 



(29) 



As in the previous section, the integration domain is the 
whole solid angle f2. Actually, when the explicit calcula- 
tions are carried out, the symmetries of the problem al- 
low restricting the 9 and tp integrals to the range [0, 7r/2]. 
The parameters that can be used for fitting the experi- 
mental results are all contained into the functions Pjj r (0) 

and P^(^) <p) appearing in the above expression 



II 

One 



of these parameters is < H c (j) >, which is the average 
of the critical fields over the weak link ensemble, and 
can be expressed in terms of the critical current intensity 
j, according to Eq. (||). Another parameter is the flux- 



compression coefficient k, which, according to Eq. (10) 



simply represents a scaling factor of the applied mag- 
netic field intensity H in the evaluation of < He{9) >. 
As already stated, the standard deviations a c and o~e can 
be taken as proportional respectively to < H c (je) > and 
< He (6) >, and the proportionality constants become 
therefore non-dimensional parameters to be used for the 
fit. Actually, as we will see later, the choice of their val- 
ues has little influence on the theoretical curves, and thus 
they can be considered as non-crucial parameters for the 
fit. As discussed in Section [v|, the most important pa- 
rameter in fitting the experimental data is < H C (J) >, 
which represents the "strength" of the weak link ensem- 
ble characteristic of the specimen under consideration. 
Results of the numerical calculations are reported and 
discussed in Section [v|. 



IV. EXPERIMENTAL RESULTS 

In this section we report the results of an extensive set 
of resistance measurements we carried out on a granular 
YBCO specimen. The sample, of about 1 x 1 x 10 mm 3 , 
was obtained by sinterization of high-purity powders 



of Y2O3, BaC03, and CuO in the stoichiometry ratio 
1:2:3. The oxygen content of the YBCO chains was 
then modified by a long-time annealing (30 days) in a 
controlled oxygen atmosphere at T = 720 K. This pro- 
cess had a twofold effect on the critical parameters of 
the material. First, it lowered the critical temperature 
down to T c o = 65 K; second, it reduced the critical cur- 
rent, that was found to be as low as 2.5T0 5 A/m 2 at 
T = 27 K and in zero field. While the first effect can be 
ascribed to a change in the oxygen doping in the grains, 
the second is mainly due to a strong weakening of the 
links between grains. This made the transport proper- 
ties of the material below T c q be mostly controlled by the 
weak links in a wide range of magnetic fields. 

The resistance measurements were carried out by us- 
ing the conventional four-probe technique. The four con- 
tacts were obtained by Ag evaporation at the opposite 
ends of the sample. Both the current and voltage leads 
were made of thin Pt wires fixed to the Ag contact by 
using Ag conductive paste. To eliminate the possible 
unwanted contributions of thermoelectric voltages, the 
current-reversal technique was used. Moreover, to avoid 
the small Joule heating of the sample, the current was 
injected into the sample only during the time strictly 
necessary for the measurement. The magnetic field was 
applied either parallel or perpendicular to the current, 
which was always flowing along the same direction, i.e. 
parallel to the longest side of the sample. 

The resistance vs. temperature curves reported in 
Fig. U were obtained at a fixed value of the magnetic 
field and different current densities (a), or at a fixed cur- 
rent and different magnetic fields (b). In both cases, the 
magnetic field H was perpendicular to the current j. 

The clear step of the curves at the temperature T c q = 
65 K indicates that, at this temperature, the grains be- 
come superconductive while all the weak links remain 
in the resistive state. It is clearly seen that the resid- 
ual resistance just below T c q is practically independent 
of both the magnetic field intensity H and the current 
density j. If the temperature is further lowered, instead, 
the R(T) curves split depending on the values of H and 
j. Fig. ^| shows the resistance vs. temperature curves 
obtained with fixed values of H and j in the two cases 
where H || j and H _L j. It is clearly seen that the resid- 
ual resistance just below T c0 is as well unaffected by the 
respective orientation of these two vectors. Since just 
below T c o practically all the weak links are still in the 
resistive state, this residual resistance can be identified 
with that of the whole weak-link network. Therefore, 
these results indirectly support our hypothesis that the 
normal-state resistance of each weak link is independent 
of H and j - and thus on their respective orientation. 

In order to study in greater detail the anisotropy of the 
resistance highlighted by the curves in Fig. ^, and to com- 
pare the predictions of our model with the experimental 
results, we measured the resistance of the sample as a 
function of the magnetic field intensity, by keeping both 
the temperature and the current fixed to a certain value. 
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FIG. 5: Experimental R vs. T curves measured in the YBCO 
specimen described in the paper, in the presence of a mag- 
netic field H perpendicular to the current density j. The 
curves in (a) were measured with a magnetic field of intensity 
H = 2.8 ■ 10 3 A/m and different values of the current den- 
sity. The curves in (b), instead, were obtained with a current 
density j = 2.8 • 10 4 A/m 2 and different values of the applied 
magnetic field. The drop of the resistance at the tempera- 
ture T c o — 65 K is due to the superconducting transition of 
the grains, while the intergrain regions remain in the resistive 
state. Notice that the residual resistance just below T c o (that 
is, the resistance of the whole weak-link network) is almost 
independent of both the current density and the magnetic 
field. 



For example, Fig. [7] shows two sets of R vs. H curves 
obtained after cooling the sample down to T = 27 K in 
zero field, and then by applying a magnetic field H per- 
pendicular to j. The two sets of curves refer to different 
values of the current density, j — 4- 10 4 A/m 2 (open and 
solid squares) and j = 6 • 10 4 A/m 2 (open and solid tri- 
angles) . For these values of the current density, self-field 
effects arc very small and thus we could neglect them. 

In Fig. 0, solid (open) symbols indicate the resistance 
measured while the magnetic field intensity is increased 
(decreased). It is clearly seen that the curves are very 
nearly reversible up to about 4.5 • 10 3 A/m. In this 
regime, the magnetic flux is very likely to be excluded 
by the grains, and the variation of the resistance due 




0.1 J 1 ■ 1 ■ 1 ' 1 ■ 1 ■ r- 

30 40 50 60 70 80 
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FIG. 6: Experimental R vs. T curves measured with j = 
2 • I0 4 A/m 2 and H = 2.8 ■ I0 3 A/m, in the two cases where 
H _L j (open triangles) and H || j (solid squares). Notice that 
the residual resistance just below T c o is independent of the 
orientation of j with respect to H. As discussed in the text, 
this suggests that the anisotropy is not due to an intrinsic 
dependence of the intergrain resistivity on the angle between 
H and j . 



to the magnetic field can be ascribed to the transition 
of weak links from the superconductive to the resistive 
state, or vice-versa. Moreover, the curves reported in 
Fig. [I] show that an approximate inverse proportionality 
exists between the magnetic field H and the current den- 
sity j for a given value of the specimen resistance. In 
other words, the same resistance, let's say R ~ 2 mfi, is 
obtained with H ~ 1.6 • 10 3 A/m and j = 6 • 10 4 A/m 2 , 
or with H ~ 2.4 • 10 3 A/m and j = 4 • 10 4 A/m 2 . This 
is exactly what is expected according to our model, if 
Eq. (^) holds true and the sample is not too close to the 
percolation threshold. 

When the field is increased above 4.5 • 10 3 A/m, the 
magnetic flux begins to penetrate into the superconduct- 
ing grains, where it remains trapped when the external 
field is decreased, thus giving rise to a strong hystere- 
sis. Let us discuss for a while what happens in this irre- 
versible regime. At the beginning of the flux penetration 
into the grains, the increment of the local magnetic field 
in the intergrain regions for a given increase of the ex- 
ternal field is smaller than in the reversible regime. This 
results in a reduction of the slope of the R versus H curve. 
For higher values of the applied field, the flux penetration 
into grains gives rise to a decrement of the local magnetic 
field, even if the external field is increased, and then the 
slope of the curve becomes slightly negative, as clearly 
shown in the figure. Within the model developed in the 
present paper, this effect can be interpreted as due to 
a reduction of the flux compression coefficient k, which 
therefore turns out to depend on the magnetic field. Ac- 
tually, we disregarded this dependence and took k as a 
constant, which is only true as long as the magnetic field 
does not penetrate into the grains. Incidentally, this is 
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FIG. 7: Resistance versus magnetic field curves at T — 27 K 
for two different values of the current density: j — 4-10 4 A/m 2 
(squares) and j = 6 • 10 4 A/m 2 (triangles). The magnetic 
field H was applied perpendicular to the current density j. 
Solid (open) symbols indicate the resistance measured when 
the field is increased (decreased), as indicated by the arrows. 
The onset of irreversibility at H m 4.5 • 10 3 A/m is due to 
flux penetration and flux trapping inside the superconduct- 
ing grains. Below this value, instead, the R vs. T curves 
are nearly reversible, which indicates the almost complete ab- 
sence of flux trapping either by grains or by superconducting 
rings. In this regime, the changes in the resistance due to the 
magnetic field are attributed to weak-link transitions. No- 
tice that the curves are in fair agreement with Eq. (^), since 
the same value of the resistance is obtained for magnetic field 
values inversely proportional to the current density j. 



one of the reasons why the validity of our model is re- 
stricted to the reversible regime. Let us just point out 
here that, if the field is further increased (as we did in 
another set of measurements not reported here) the slope 
becomes positive again. At about 1.6 • 10 5 A/m, the re- 
sistance becomes practically constant and saturates to 
a current-independent value. Notice that the same be- 
haviour can be observed, close to T c q, in the curves re- 
ported in Fig. ||a. 

In Fig. U two R versus H curves are shown, obtained 
with the same values of the current density (4- 10 4 A/m 2 ) 
but in the two cases H _L j and H || j. By starting from 
these data sets, the magnetic field- dependence of the ra- 
tio i] — Ru/Rx can be easily obtained. As a matter of 
fact, this dependence is reported in Fig. |9| (solid trian- 
gles) together with a similar curve obtained with a cur- 
rent density j = 6 • 10 4 A/m 2 (solid squares). In the 
same figure, the best-fitting curves calculated by using 
Eq. (E9) are also shown (open circles). The numerical 
integration of Eq. (^9|) was performed by means of the 
computer program Macsyma V. 2.2, by Macsyma Inc. 
The values of the best-fit parameters, reported in the fig- 
ure caption, will be discussed in the next section. Let 
us just point out here that, in spite of the many simpli- 
fications implicit in our model, there is a fair agreement 



FIG. 8: Resistance vs. magnetic field curves at the same 
values of current density and magnetic field, in the two cases 
where H _L j (triangles) and H || j (circles). The anisotropy 
of the resistance induced by the magnetic field is clearly seen. 
The values of 77 reported in Fig. ^ are taken from these curves 
and from a similar couple of curves measured with a current 
density j = 6 • 10 4 A/m 2 . 

between theoretical and experimental data. 

V. DISCUSSION AND CONCLUSIONS 

The results reported in Fig. [)] show that there is a gen- 
eral agreement between the experimental data and the 
results of our model. The agreement is particularly good 
for intermediate values of the magnetic field intensity. 
Instead, some deviations between experimental and the- 
oretical results are evident for low and high values of the 
applied field. On account of what was already pointed 
out in the previous sections, these deviations are not sur- 
prising. As a matter of fact, our model has been devel- 
oped under the assumption that the specimen was neither 
too close to the superconductive percolation threshold 
(that takes place when the magnetic field is too weak), 
nor in the condition of flux penetration within the grains 
(that occurs when the field is too intense). 

Let us now focus on the best-fitting values of the pa- 
rameters that enter in our model. In principle, the ad- 
justable parameters of the model are: the mean critical 
field < H c (j) >, the flux compression factor k, and the 
standard deviations og and a c , characterizing the spread 
of the local values of Hi(9) and H c (j). Since, according 
to the assumptions made above, a a and a c are taken as 
being proportional to < Hg(9) > and < H c (j) > respec- 
tively, it can be easily shown that P^ T (0) and P±(0,ip) 
depend mainly on the quantity H*(j) =< H c (j) > /k 
and, to a minor extent, on the proportionality factors for 
a i and a c . H*{j) represents the reduced mean critical 
field of the weak links, which depends on the intensity of 
the current density, j, and on the flux compression fac- 
tor, k. Thus, according to the present model, comparing 
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FIG. 9: Comparison between experimental and theoretical 
results concerning the anisotropy factor r\ for two different 
values of the macroscopic current density j. Solid symbols 
represent experimental data, while open circles are theoretical 
points calculated by means of Eq.(E9). The theoretical curves 
shown here are those that best fit the experimental data, and 
were obtained by taking H*(j) = 3.0 • 10 3 A/m and H*(j) = 
2.0- 10 3 A/m for j = 4-10 4 and 6-10 4 A/m 2 , respectively. The 
ratio between the two values of H*(j) is clearly the inverse 
of the ratio between the current densities, in agreement with 
Eq. (^). The widths of the distributions of critical fields and 
local fields were taken as being a c — 0.3- < H c (je) > and 
ae = 0.7- < He(9) > for both the curves. Notice that the 
actual values of < H c (je) > may be an order of magnitude 
greater than those of H*(j), owing to the effect of the flux 
compression coefficient k. 

the experimental data with the theoretical results allows 
the determination of H*(j) as the main best- fit param- 
eter, whose value is only little affected by the other two 
parameters ae and a c . H c (j) can thus be taken as a 
quantity apt to characterize the "strength" of the weak- 
links ensemble and it is strongly dependent on the type of 
granular superconducting material. It has been observed 
that thermal treatments can produce strong changes in 
this quantity, leaving practically unaffected the super- 
conducting properties of the grains, i.e. the critical tem- 
perature and the value of the magnetic field at which the 



flux begins to penetrate into the grain; 




In conclusion, we have presented a simple model which 
is able to explain the anisotropy of the resistance shown 
by granular HTSCs in the presence of a magnetic field. 
The model can be applied if the intensity of the applied 
magnetic field is such that the bulk superconductivity is 
disrupted, but the grains remain in the superconducting 
state. This physical requirement is easily fulfilled in sam- 
ples where the connections between grains is weak, as a 
result, for instance, of thermal treatments and anneal- 
ing processes. In these conditions, we have shown that 
there is no need of assuming an intrinsic dependence of 
the intergrain resistivity on the angle between H and j, 
as instead the models based on the conventional theories 
of current-driven flux motion do. Incidentally, the mo- 
tion of unpinned vortices is actually the main origin of 
the resistance in those materials where the links between 
grains are so strong that the flux can penetrate into the 
grains and even become unpinned well before the resistive 
transition of the weak links sets irfcl. 

If the material is made up of weakly linked grains, in- 
stead, the local magnetic field and the local current den- 
sity in a given weak link are always perpendicular to each 
other, provided that the flux penetration into the grains 
is negligible. In this case the proposed model, which de- 
scribes the anisotropy of the resistivity as being due to 
the spatially anisotropic distribution of the resistive weak 
links seems more appropriate. Actually, the results of the 
model are in reasonable agreement with the experimental 
data obtained from resistance measurements in a YBCO 
granular specimen which satisfies the above mentioned 
conditions. Good experimental evidence is also found 
for the inverse proportionality between the mean critical 
field of the weak links and the current density, which is 
the basic assumption over which the proposed model has 
been developed. 
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